Phase-Space analysis of Teleparallel Dark Energy 
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We perform a detailed dynamical analysis of the teleparallel dark energy scenario, which is based 
on the teleparallel equivalent of General Relativity, in which one adds a canonical scalar field, 
allowing also for a nonminimal coupling with gravity. We find that the universe can result in the 
quintessence-like, dark-energy-dominated solution, or to the stiff dark-energy late-time attractor, 
similarly to standard quintessence. However, teleparallel dark energy possesses an additional late- 
time solution, in which dark energy behaves like a cosmological constant, independently of the 
specific values of the model parameters. Finally, during the evolution the dark energy equation- 
of-state parameter can be either above or below —1, offering a good description for its observed 
dynamical behavior and its stabilization close to the cosmological-constant value. 
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I. INTRODUCTION. 
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According to observations, the observable universe ex- 
periences an accelerating expansion at late times 
There are two major ways to explain such a behavior, 
apart from the simple consideration of a cosmological 
constant. The first is to modify the gravitational sec- 
tor itself (see for a review and references therein), 
obtaining a modified cosmological dynamics. The other 
approach is to modify the content of the universe intro- 
ducing the dark energy with negative pressure, which can 
be based on a canonical scalar field (quintessence) [1, 0] , 
a phantom field [H| , or on the combination of both these 
fields in a unified scenario called quintom Q (see the 
reviews and 3). 

"Teleparallel" dark energy is a recently proposed sce- 
nario that tries to incorporate the dark energy sector 

[l^. It is based on the "teleparallel" equivalent of 
General Relativity (TEGR), that is based on torsion in- 
stead of curvature formulation [111, [i3 1 i^. which one adds 
a canonical scalar field, and the dark energy sector is 
attributed to this field. In the case where the field is 
minimally-coupled to gravity, the scenario is completely 
equivalent to standard quintessence 0] , both at the back- 
ground and perturbation levels. However, in the case 
where one switches on the nonminimal coupling between 
the field and the torsion scalar which is the only suitable 
gravitational scalar in TEGR, the resulting scenario has a 
richer structure, exhibiting quintessence- like or phantom- 
like behavior, or experiencing the phantom-divide cross- 
ing 0. 

Since teleparallel dark energy exhibits interesting cos- 
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mological behavior, in the present work we perform a 
phase-space and stability analysis of such a scenario, ex- 
amining in a systematic way the possible cosmological be- 
haviors, focusing on the late-time stable solutions. Such 
an approach allows us to bypass the high non-linearities 
of the cosmological equations, which prevent any com- 
plete analytical treatment, obtaining a (qualitative) de- 
scription of the global dynamics of these models, that 
is independent of the initial conditions and the specific 
evolution of the universe. Furthermore, in these asymp- 
totic solutions we calculate various observable quantities, 
such as the deceleration parameter, the effective (total) 
equation-of-state parameter, and the various density pa- 
rameters. 

The plan of the work is the following: In section HIl we 
briefly review the scenario of teleparallel dark energy and 
in section mil we perform a detailed phase-space analysis. 
In section IIVI we discuss the cosmological implications 
and the physical behavior of the scenario. Finally, section 
fVlis devoted to the summary of the results. 



II. TELEPARALLEL DARK ENERGY 

In this section we review teleparallel dark energy. Such 
a scenario is based on the "teleparallel" equivalent of 
General Relativity (TEGR) [ill. El, in which instead of 
using the torsionless Levi-Civita connection one uses the 
curvatureless Wcitzenbock one. The dynamical objects 
are the four linearly independent, parallel vector fields 
called vierbeins, and the advantage of this framework is 
that the torsion tensor is formed solely from products of 
first derivatives of the tetrad. In particular, the vierbein 
field eA{x^) forms an orthonormal basis for the tangent 
space at each point x'^, that is ba ■ gb = i]ab^ where 
rjAB = diag{\,—\,—\,—l), and furthermore the vector 
ba can be analyzed with the use of its components 
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in a coordinate basis, that is 



=19 



In such a 



construction, the metric tensor is obtained from the dual 
vierbein as 



= VAB e^{x) ef (x). 



(1) 



Furthermore, the torsion tensor of the Wcitzenbock con- 
nection r^,^ [l^ reads 



rj-i/' 



p, J- 



(2) 



In the present formahsm all the information concern- 
ing the gravitational field is included in the torsion tensor 
T^^ . The corresponding "teleparallel Lagrangian" can be 
constructed from this torsion tensor under the assump- 
tions of invariance under general coordinate transforma- 
tions, global Lorentz transformations, and the parity op- 
eration, along with requiring the Lagr angi an density to 
be second order in the torsion tensor jl2| . In particu- 
. ar. it is proved to be just the torsion scalar T, namely 

ni[ii[ii: 

C^T= ir'""^Tp,„. + ^T'^-'n^, - Tp/T"^, . (3) 

Therefore, the simplest action in a universe governed by 
teleparallel gravity is 



d xe 



T 

2^ 



(4) 



where e = det(e^) = \/—g (one could also include a 
cosmological constant). Variation with respect to the 
vierbein fields gives the equations of motion, which are 
exactly the same as those of General Relativity for every 
geometry choice, and that is why the theory is called 
"teleparallel equivalent to General Relativity" . 

In principle one has two ways of generalizing the action 
(ID), inspired by the corresponding procedures of stan- 
dard General Relativity. The first is to replace T by an 
arbitrary function /(T) [l^[i^, similarly to f{R) exten- 
sions of GR, and obtain new interesting terms in the field 
equations. The other, on which we focus in this work, is 
to add a canonical scalar field in similarly to GR 
quintessence, allowing for a nonminimal coupling with 
gravity. Then the dark energy sector will be attributed 
to this field and therefore the corresponding scenario is 
called "teleparallel dark energy" Q. In particular, the 
action will simply read: 



S 



T 

2^ 



(5) 



^ We follow the notation of [T^, that is Greek indices fj.,u,... and 
capital Latin indices A,B,... run over all coordinate and tan- 
gent space-time 0, 1, 2, 3, while lower case Latin indices (from 
the middle of the alphabet) ... and lower case Latin indices 
(from the beginning of the alphabet) a, b,... run over spatial and 
tangent space coordinates 1, 2, 3, respectively. Finally, we use 
the metric signature (-|-, — , — , — ). 



Concerning the nonminimal coupling we emphasize that 
since in TEGR, that is in the torsion formulation of GR, 
the only scalar is the torsion one, the nonminimal cou- 
pling will be between this and the scalar field (similarly 
to standard nonminimal quintessence where the scalar 
field couples to the Ricci scalar). 

Variation of action ([5]) with respect to the vierbein 
fields yields equation of motion 
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P rp V 



(6) 



where T stands for the usual energy-momentum ten- 
sor. Therefore, imposing a flat Friedmann-Robertson- 
Walker (FRW) background metric 



ds^ ^ dt'^ - a^{t) Sijdx^'dx^ , 



(7) 



where t is the cosmic time, a;* are the comoving spatial 
coordinates and a{t) is the scale factor, that is for a vier- 
bein choice of the form 



= diag(l,a,a,a), 
we obtain the corresponding Friedmann equations: 
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H 



K 



P0 + Pv 



P<j,+Pcli+ Prn + Pn 



(8) 

(9) 
(10) 



where H = a/a is the Hubble parameter, a dot denotes 
differentiation with respect to t. In these expressions, 
Pm and Pm are the matter energy density and pressure, 
respectively, following the standard evolution equation 
Pm + 3-ff(l + Wm)Pm = 0, with Wrn = Pm/Pm the mat- 
ter cquation-of-state parameter. Furthermore, we have 
introduced the effective energy density and pressure of 
scalar field 
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(11) 



V{(f>) + A^H(f>,p + n3H^ + 2H] (12) 
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Moreover, variation of the action with respect to the 
scalar field provides its evolution equation, namely: 



4> + + 6^H^(j) + l/'(0) = 0. 



(13) 



Note that in the above expressions we have used the use- 
ful relation 



6i^^ 



(14) 



which straightforwardly arises from the calculation of ([3]) 
for the fiat FRW geometry. 
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In this scenario, similarly to the standard quintessence, 
dark energy is attributed to the scalar field, and thus its 
equation-of-state parameter (wde) reads: 



P4> 

Wde = w^ = — . 



(15) 



Finally, one can see that the scalar field evolution 
leads to the standard relation 



/50 + 3H{1 + w^)p^ = 0. 



(16) 



Teleparallel dark energy exhibits very interesting cos- 
mological behavior. In the minimally-coupled case the 
cosmological equations coincide with those of the stan- 
dard quintessence, both at the background and pertur- 
bation levels. However, in the nonminimal case one 
can obtain a dark energy sector being quintessence-like, 
phantom-like, or experiencing the phantom-divide cross- 
ing during evolution, a behavior that is much richer com- 
paring to General Relativity (GR) with a scalar field 
[3, ll^. Therefore, it is both interesting and necessary 
to perform a phase-space analysis, that is to investigate 
late-time solutions that are independent from the initial 
conditions and the specific universe evolution. This is 
performed in the next section. 



III. PHASE-SPACE ANALYSIS 



In a. Then we extract its critical points Xc satisfying 
X' = 0, and in order to determine the stability proper- 
ties of these critical points, we expand around Xc, setting 
X = Xc -l- U, with U the perturbations of the variables 
considered as a column vector. Thus, for each critical 
point we expand the equations for the perturbations up 
to the first order as: U' = Q • U, where the matrix Q 
contains the coefficients of the perturbation equations. 
Finally, for each critical point, the eigenvalues of Q de- 
termine its type and stability. 

In the scenario at hand, we introduce the auxiliary 
variables: 



y 



(17) 



Using these variables the Friedmann equation © be- 
comes 



x2+y2-z2sgn(0 + ^ = l. 



(18) 

This constraint allows for expressing p,„ as a function of 
the auxiliary variables ((T7| . Therefore, using (fT8)) and 
(1111) we can write the density parameters as: 



In order to perform the phase-space and stability anal- 
ysis of the scenario at hand, we have to transform the 
aforementioned dynamical system into its autonomous 
form X' = f(X) [H [Tsj . where X is the column vec- 
tor constituted by suitable auxiliary variables, f(X) the 
corresponding column vector of the autonomous equa- 
tions, and prime denotes derivative with respect to M = 



I - X - y + z sgn(^) 



nuE^^=x^+v'-z\gn{£,), (19) 

while for the dark-energy equation-of-state parameter 
([T5|) we obtain: 



2 _ „,2 , /I /2 /Tr\a„^(C\ _ ^2 



[1 + z2 sgn(C)] + y2 _ ^2sgn(^)] ■ (20) 



In this scenario wde can be quintessence-like, phantom- 
like, or experience the phantom divide crossing during 
the evolution. Without loss of generality, in the following 
we restrict the analysis in the dust matter case, that is 
we assume that Wm = 0. 

It is convenient to introduce two additional quanti- 
ties with great physical significance, namely the "total" 
equation-of-state parameter: 



P<t> 



P4, + Pv 



wdeO.de 



y'+4y'fz2:\/le|sgn(e) 
1-1-2:2 sgn(^) 

(21) 



and the deceleration parameter q: 



H 1 3 

-1 = 1 Wfnt 

i/2 2 2 



1 + 3(^2 -2/2)+ z-|-4y6a;VW 2 sgn(0 



-.(22) 



2 [1 + z2 sgn(0] 

Finally, concerning the scalar potential V{4)) the usual 
assumption in the literature is to assume an exponential 
potential of the form 



V = Voexp(— kA0), 



(23) 
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since exponential potentials are known to be significant 
in various cosmological models [13, [13 (equivalcntly, but 
more generally, we could consider potentials satisfying 
~ ~ KV{<t>) ^00^'' ^ const, which is valid for arbitrary 
but nearly flat potentials [l^). Moreover, note that the 
exponential potential was used as an example in the ini- 
tial work on teleparallel dark energy [^ . 

In summary, using the auxiliary variables (|17p and con- 
sidering the exponential potential (|23p . the equations of 
motion ([9]), ([T0|) and (fT3|) in the case of dust matter, can 
be transformed into the following autonomous system: 



2z2sgn(e) + 2 
1 



2 



2^ 



Vex 



z2sgn(0 + 2 
3x 



sgn(^)z2 + I 



2sgn(C)z2 + 2 2 



1 - 



z2sgn(^) + 1 







zsgn(^) + a| 


1 + z'^sgn{S,) 



(24) 



where we have used that for every quantity F we acquire 
F = HF' . Since Pm is nonncgative, from ((T5)) we obtain 
that + y'^ — z^sgn(^) < 1. Thus, we deduce that for 
^ < the system ([M]) defines a flow on the compact phase 
space ^ = {a;^ + — z^sgn(^) < 1, ?/ > O} , for ^ = 
the phase space is compact and it is reduced to the circle 
\[/ = -I- 2/^ < 1, J/ > O} , while for ^ > the phase 
space ^ is unbounded. 

Before proceeding we make two comments on the de- 
grees of freedom and the choice of auxiliary variables. 
Firstly, as we have already mentioned, in the minimal 
coupling case, that is when ^ = 0, the model at hand co- 
incides with standard quintessence, which phase space 
analysis is well known using two degrees of freedom, 
namely the variables x and y defined above [l^- On 
the other hand, when ^ ^ we have three degrees of 
freedom and all a;, y, z are necessary. Therefore, in order 
to perform the analysis in a unified way, we use the three 
variables defined above, having in mind that for ^ = the 
variable z becomes zero and thus irrelevant. Secondly, 
apart from the standard choices of the variables x and y, 
one must be careful in suitably choosing the variable z 
in order not to lose dynamical information. For example, 
although the choice for the variable z = is another 
reasonable choice, however it still loses the critical points 
that lie at "infinity". Therefore, in order to completely 
cover the phase-space behavior in the following we will 
additionally use the Poincare central projection method 
to investigate the dynamics at "infinity" . The negligence 
of this point was the reason of the incomplete phase space 
analysis of teleparallel dark energy performed in [20| . 



A. Finite Phase-space analysis 

Let us now proceed to the phase-space analysis. The 
real and physically meaningful (that is corresponding to 
an expanding universe, i.c possessing H > 0) critical 
points {xc,yc,Zc) of the autonomous system pi]) , ob- 
tained by setting the left hand sides of the equations to 
zero, are presented in Table HI In the same table we pro- 
vide the conditions for their existence. The 3x3 matrix 
Q of the linearized perturbation equations of the system 
([Ml) is shown in Appendix |^ For each critical point of 
Table U we examine the sign of the real part of the eigen- 
values of Q in order to determine the type and stability 
of the point. The details of the analysis and the various 
eigenvalues are presented in Appendix [Xj and in Table |T] 
we summarize the results. ^ In addition, for each critical 
point we calculate the values of floE, wde, Wtot and q 
given by ^ and ([22]) respectively. 



B. Phase-space analysis at infinity 



Due to the fact that the dynamical system ([24)) is non- 
compact for the choice ^ > 0, there could be features 
in the asymptotic regime which are non-trivial for the 
global dynamics. Thus, in order to complete the analysis 
of the phase space we must extend our study using the 
Poincare central projection method [2l|. 

Let us consider the Poincare variables 



Xr ~ /9 COS 6* sin -0, z,. = p sin 6* sin?/', y,. = pcosij}, (25) 

7/2 + e e [0,27r], and 



where p ~ 



r 



,,2 



~f < ^ f (we restrict the angle ip to this interval 
since the physical region is given by y > 0). Thus, the 
points at "infinity" (r — >■ -|-oo) are those having p ^ 1. 
Furthermore, the physical phase space is given by 



{Xr,yr,Zr) G [-1,1] X [0,1] X [-1,1] < + < 1 



Inverting relations ([25]) and substituting into ([19)) . ([20]) . 
we obtain the dark energy density and cquation-of- 
statc parameters as a function of the Poincare variables, 
namely: 



^DE — 



WDE 



(3x^-|-4^/6^/^2^Xr-3y^)(l- 



3(1 — x^. — ?/p)(a:^-}-y2_^2j 



(26) 
(27) 



^ The critical point D that exists only for ^ = 0, is very close 
to be a stationary point in the general case 0) too, but 

with an arbitrarily small value of the slope of the potential A. It 
becomes an exact critical point only for A = 0, namely the point 
J. Furthermore, for arbitrary large values of A the singular point 
E (that exists for § = 0) mimics the behavior of the point A (that 
exists for ^ ^ 0). 
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Cr. P. 


Xc 


Vc 


2;c 


Exists for 


Stability 


^DE 


WDE 


Wtot 


q 


A 











all C,A 


saddle point 





arbitrary 





i 


B 


1 








^ = 0, all A 


unstable for A < \/6 
saddle point otherwise 


1 


1 


1 


2 


C 


-1 








^ = 0, all A 


unstable for A > — \/6 
saddle point otherwise 


1 


1 


1 


2 


D 


A 

V6 







C = 0, A^ < 6 


stable node for A^ < 3 
saddle point for 3 < A^ < 6 


1 






-1 + ^ 

^2 


E 


7^ 


7^ 

V Ix 





C = 0, A^ > 3 


stable node for 3 < A^ < ^ 
stable spiral for A^ > ^ 


3 








1 
2 


F 







AC 


^/i€^ 


< A^ < ^ 


stable node for A^ < ^ 


1 


-1 


-1 


-1 


G 





y2€+2\/«(5-A2) 


5+\/«(S-^") 
AC 




0<A^<^or$<0 


saddle point 


1 


-1 


-1 


-1 


J 





1 





^ / and A = 


stable spiral for | < C 
stable node for < ^ < | 
saddle point for ^ < 


1 


-1 


-1 


-1 



TABLE I: The real and physically meaningful critical points of the autonomous system (|24|l , and the corresponding values of 
the dark-energy density parameter Q.de, of the dark-energy equation-of-state parameter wde, of the total equation-of-state 
parameter wtot and of the deceleration parameter q. 



Cr. P. 


Xrc 


yrc 




Stability 


^Ide 


WDE 


Wtot 


q 


K± 


+ 2 





1 \/2 
^ 2 


unstable for < C < | 
saddle point for | < C 


arbitrary 


arbitrary 




2-2\/6C 


L± 


1 V2 

^ 2 





<V2 

^ 2 


saddle point for all ^ > 


arbitrary 


arbitrary 




2 + 2^ 



TABLE IL The real and physical critical points of the autonomous system (|24p at infinity, which exists for ^ > 0, all A, and 
the corresponding values of the dark-energy density parameter ^Ide, of the dark-energy equation-of-state parameter wde, of 
the total equation-of-state parameter wtot and of the deceleration parameter q. 



and similarly substituting into ([2T|). p2|) we obtain the 
corresponding expressions for the total equation-of-state 
and deceleration parameters: 



Wtot 



3a;2 + 4^67^, 



3(1 



(28) 



correspond to the limit 

\(\)\ oo, 

satisfying the rate (ln(/))' = \f^- = a/6 cot ( 



(30) 



Applying the procedure prescribed in Appendix IB] we 
conclude that there are four physical critical points at in- 
finity, namely the points K±^ satisfying cot^ = — 1, and 
the points L± satisfying cotfl = 1. These critical points, 
along with their stability conditions are presented in Ta- 
ble ini In the same Table wc include the corresponding 
values of the obscrvablcs fi_D_E, wde^ Wtot and q, calcu- 
lated using (|26l).([27|.p8|.(f29l). 

As we show in Appendix |BJ the above critical points 



IV. COSMOLOGICAL IMPLICATIONS 

Having performed the complete phase-space analysis 
of teleparallcl dark energy, we can now discuss the cor- 
responding cosmological behavior. A first remark is that 
in the minimal case (that is ^ = 0) we do verify that 
the scenario at hand coincides completely with standard 
quintessence. Therefore, we will make a brief review on 
the subject and then focus on the nonminimal case. 

The points B to E exists only for the minimal case, 
that is only for ^ = 0. Points B and C arc not stable, 
corresponding to a non-accelerating, dark-energy domi- 
nated universe, with a stiff dark-energy equation-of-state 
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parameter e qual to 1. Both of them exist in standard 
quintessence [l7l |. 

Point I? is a saddle one for 3 < < 6, however for 
< < 3 it is a stable node, and thus it can attract the 
universe at late times. It corresponds to a dark-energy 
dominated universe, with a dark-energy equation-of-state 
parameter lying in the quintessence regime, which can 
be accelerating or not according to the A- value. This 
point exists in standard quintessence |17| . It is the most 
important one in that scenario, since it is both stable and 
possesses a wde compatible with observations. 

Point E is a stable one for 3 < A^. It has the ad- 
vantage that the dark-energy density parameter lies in 
the interval < ftoE < ^, that is it can alleviate the 
coincidence problem, since dark matter and dark energy 
density parameters can be of the same order (in order 
to treat the coincidence problem one must explain why 
the present dark energy and matter density parameters 
are of the same order of magnitude although they fol- 
low different evolution behaviors). However, it has the 
disadvantage that wde is and the expansion is not ac- 
celerating, which are not favored by observations. This 
point exists in standard quintessence [T7| . 

Let us now analyze the case ^ 7^ 0. In this case we 
obtain the critical points A, F, G and J. The point A is 
saddle point, and thus it cannot be late-time solution of 
the universe. It corresponds to a non-accelerating, dark- 
matter dominated universe, with arbitrary dark-energy 
equation-of-state parameter. Note that this trivial point 
exists in the standard quintessence model too [17[ , since 
it is independent of ^. 

The present scenario of teleparallel dark energy, pos- 
sesses two additional, non-trivial critical points that do 
not exist in standard quintessence. Thus, they account 
for the new information of this richer scenario, and as 
expected they depend on the nonminimal coupling ^. In 
particular, point F is stable if A^ < ^, and thus it can 
attract the universe at late times. It corresponds to an 
accelerating universe with complete dark energy domina- 
tion, with Wde = —1, that is dark energy behaves like 
a cosmological constant. We stress that this wde value 
is independent of A and ^, which is an important and 
novel result. Thus, while point D (the important point 
of standard quintessence) needs to have a very flat, that 
is quite tuned, potential in order to possess a wde near 
the observed value —1, point F exhibits this behavior for 
every A-value, provided that A^ < ^. This feature is a sig- 
nificant advantage of the scenario at hand, amplifying its 
generality, and offers a mechanism for the stabilization of 
Wde close to the cosmological-constant value. Similarly, 
we have the point G, which has the same cosmological 
properties with F, however it is a saddle point and thus 
it cannot be a late-time solution of the universe, but the 
universe can spend a large period of time near this solu- 
tion. 

Finally, when £, and for the limiting case A = 0, 
that is for a constant potential, the present scenario ex- 
hibits the critical point J, which is stable for ^ > 0. It 



corresponds to a dark-energy dominated, de Sitter uni- 
verse, in which dark energy behaves like a cosmological 
constant. 

Let us now now analyze the critical points at infinity. 
The points K± are saddle for ^ > |. The fact that they 
possess arbitrary VIde and wde is a significant advan- 
tage, since such solutions can alleviate the coincidence 
problem. Note however that the corresponding Wtot and 
deceleration parameter arc not arbitrary, offering a good 
quantitative description of the cosmological behavior (ac- 
tually this is the reason we introduced these observables). 
In particular, we conclude that for ^ > ^ we obtain an 
accelerated expansion, while for C > f the expansion is 
super-accelerating (q < — 1 that is H > and wtot < ^1) 
that is the universe presents a phantom behavior. We 
mention that since these points arc saddle they cannot 
attract the universe at late times, however the universe 
can spend a large period of time near these solutions, 
before approaching the saddle point A or the global at- 
tractor F (for potentials with slope A^ < 0- that case 
the above physical features present a transient character, 
which is quite significant from the observational point of 
view. 

Points L± are saddle for all ^ > 0. They have arbi- 
trary Qde and wde however they correspond to a non- 
accelerating expansion, and thus they could be important 
only as transient states of the universe. 

In order to present the aforementioned behavior more 
transparently, we first evolve the autonomous system ([M]) 
numerically for the choice A = 1, and ^ = 0, that is 
in the minimal scenario, and in Fig. [T] we depict the 
corresponding phase-space behavior, projected in the x — 
y plane. As we can wee, in this case the quintessence-like 
critical point D is the late-time solution of the universe. 
Note that this point exists only for the minimal case (^ = 
0), where our model coincides with standard quintessence 
(see Fig. 2 of [13). 

As we mentioned above, the late-time attractors D and 
E exists only for the case ^ = 0. To illustrate how sensi- 
ble is the scenario at hand in slight changes of ^, in Fig. [5] 
we depict the projection of the phase-space evolution on 
the x — y plane, for A = 2 and ^ = —10"'^. In this case the 
point E becomes "quasi-stationary" , since typical trajec- 
tories remain very close to it for a large but finite time 
interval, while D becomes also quasi-stationary, but it is 
unstable to perturbations in the z-axis. ^ 

Finally, to illustrate the dynamics for positive ^, which 
is expressed through points at "infinity" . in Fig. [3] 
we present the three-dimensional Poncarc (global) phase 
space for A = 0.5 and ^ = 1. According to Table lU the 
critical points K± are saddle ones and exhibit phantom 
behavior, while L± are saddle with a stiff dark energy sec- 



^ Although the singular points B, C, D and E do not exist as criti- 
cal points for the choices ^ ^ 0, in Fig. [2]wc mark their locations 
explicitly, for comparison with the case g = of Fig. [T] 
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FIG. 1: The projection of the phase-space evolution on the 
x—y plane, for the teleparallel dark energy scenario with A = 1 
and 5 = 0. The trajectories are attracted by the quintessence- 
like stable point D, which exists only in the minimal case — 
0). 




0,0 



FIG. 3: Poncare (global) phase space for the teleparallel dark 
energy scenario with A = 0.5 and ^ = 1. The points at infinity 
K± , L± are saddles. The trajectories on both finite and infi- 
nite regions are finally attracted by the cosmological-constant- 
like stable point F (global attractor), having passed through 
matter- dominated solutions or solutions with comparable dark 
matter and dark energy sectors. 



In 




X 



FIG. 2: The projection of the phase-space evolution on the 
x—y plane, for the teleparallel dark energy scenario with A = 2 
and 5 = —10^''. Observe that an open set of trajectories 
remains very close to the quasi-stationary solution E, before 
approaching an invariant arc above the saddle point A. The 
critical point G is a saddle one. 



tor. Thus, far from their basm of attraction (but inside 
the invariant set yr = 0) the orbits departs from them 
to approach the matter-dominated solution A, while for 
Ur > Q the orbits are attracted by the cosmological- 
constant-like solution F. Thus, the epoch sequence 
K_ L_ F represents the transi- 

tion from a universe with a phantom-like dark energy, to 
a matter-dominated universe with non-phantom dark en- 
ergy, and then to a dark-energy-dominated, cosmological- 
constant-like solution. This sequence has a great cos- 
mological significance, since it can describe the epoch 
sequence inflation, dark-matter domination, dark-energy 
domination, in agreement with observations. Addition- 
ally, the epoch sequence L_ — > A'+ F represents the 
transition from a universe with comparable dark matter 
and dark energy sectors, to a cosmological-constant-like 
solution. 

Before closing this section, let us make a comment 
on another crucial difference of teleparallel dark energy, 
comparing with standard quintessence, that is of the 
^ 7^ case comparing to the ^ = one. In particular, 
when ^ = 0, in which teleparallel dark energy coincides 
with standard quintessence, wde is always larger that 
— 1, not only at the critical points, but also through- 
out the cosmological evolution as well. However, for 
^ 7^ 0, during the cosmological evolution wde can be 
either above or below —1, and only at the stable criti- 
cal point it becomes equal to —1. Such a cosmological 
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behavior is much richer, and very interesting, both from 
the theoretical as well as from the observational point of 
view, since it can explain the dynamical behavior of wde 
either above or below the phantom divide, and moreover 
its stabilization to the cosmological-constant value. 








FIG. 4: The phase-space evolution for the teleparallel dark 
energy scenario with A = 0.7 and ^ = 1. The trajecto- 
ries are attracted by the cosmological-constant-like stable point 
F, however the evolution towards it possesses a wde being 
quintessence-like, phantom-like, or experiencing the phantom- 
divide crossing, depending on the specific initial conditions. 
The orbits with a thick gray curve on the top left are those 
with Wde > —1 initially, while the thin black curve corre- 
sponds to Wde < — 1 initially. 

In order to present the novel features of the scenario at 
hand, in Fig. 2] we depict the three-dimensional phase- 
space behavior for the choice A = 0.7 and ^ = 1. In 
this case the universe at late times is attracted by the 
cosmological-constant-like stable point F. However, dur- 
ing the evolution towards this point the dark-energy 
equation-of-state parameter wde presents a very inter- 
esting behavior, and in particular, depending on the 
specific initial conditions, it can be quintessence-like, 
phantom-like, or experience the phantom-divide cross- 
ing. Such a behavior is much richer than standard 
quintessence, and reveals the capabilities of teleparallel 
dark energy scenario. 

V. CONCLUSIONS 

In the present work we investigated the dynamical be- 
havior of the recently proposed scenario of teleparallel 



dark energy [§, which is based on the teleparallel 
equivalent of General Relativity (TEGR)jthat is on its 
torsion instead of curvature formulation fill , [l^ . In this 
model one adds a canonical scalar field, in which the dark 
energy sector is attributed, allowing also for a nonmin- 
imal coupling between the field and the torsion scalar. 
Thus, although the minimal case is completely equiv- 
alent with standard quintessence, the nonminimal sce- 
nario has a richer structure, exhibiting quintessence-like 
or phantom-like behavior, or experiencing the phantom- 
divide crossing [l^ . 

Performing a detailed phase-space analysis of telepar- 
allel dark energy, we saw that the standard quintessence- 
like, late-time solution exists only for ^ = 0, and 
we showed that this solution becomes quasi-stationary 
if ^ is perturbed even slightly. The same results hold 
for the other standard quintessence solution, namely the 
stiff dark-energy late-time attractor, in which dark en- 
ergy and dark matter can be of the same order, thus 
alleviating the coincidence problem. 

Apart from the above late time solutions that exist 
also in the standard quintessence scenario, teleparallel 
dark energy has an additional and physically very inter- 
esting late-time behavior. In particular, it possesses a 
late-time attractor in which dark energy behaves like a 
cosmological constant, independently of the specific val- 
ues of the model parameters (provided that the nonmin- 
imal coupling ^ is larger than the square of the poten- 
tial exponent). This feature is a significant advantage of 
the scenario at hand, amplifying its generality, since it 
provides a natural way for the stabilization of the dark 
energy equation-of-state parameter to the cosmological 
constant value, without the need for parameter-tuning. 

Additionally, we showed that teleparallel dark energy 
admits critical points at "infinity" , which are saddle for 
^ > |, and thus although they cannot be the late-time 
solutions, the universe can spend a large period of time 
near them. Thus, one can obtain the transition from a 
universe with a phantom-like dark energy, to a matter- 
dominated universe with non-phantom dark energy, and 
then to a dark-energy-dominated, cosmological-constant- 
like solution. This sequence has a great cosmological sig- 
nificance, since it can describe the epoch sequence infla- 
tion, dark-matter domination, dark-energy domination, 
in agreement with observations 

Finally, perhaps the most interesting feature of telepar- 
allel dark energy is that in the nonminimal case the dark 
energy equation-of-state parameter can be either above 
or below —1, and only at the stable critical points it 
becomes equal to —1. Such a cosmological behavior is 
much richer comparing to standard quintessence, and 
very interesting, both from the theoretical as well as from 
the observational point of view, since it can explain the 
dynamical behavior of wde either above or below the 
phantom divide, and moreover its stabilization to the 
cosmological-constant value. All these features make the 
teleparallel dark energy a good candidate for the descrip- 
tion of dark energy. 
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Appendix A: Stability of the finite critical points 

For the critical points {xc,yc,Zc) of the autonomous 
system the coefficients of the perturbation equa- 

tions form a 3 X 3 matrix Q, which reads: 



Qii 

Ql2 
Ql3 

Q21 
Q22 

Q23 

Q31 
Q32 
Q33 



9a;2 - 3y2 4^/6.Tzy^^sgn(0 3 
2 + 2z2sgn(0 ^ 2 + z2sgn(e) 2 



\/6A- 



3a; 



1 + z2sgn(^) 



2V6 [z^sgn(0 - 2] y^x^sgn(g) ^ i {y^ - x^) zx sgnjQ _ /g 

[2 + z2sgn(e)]2 [l + zHgn{£,)f Vl^lsgn^J 



ixy 



Az2sgn(e) - 4z Vie|sgn(0 + A 



1 4- z2sgn(^) 1 + z2sgn(^) 

3 + 3a-2 - V6Aa- - 92/2 + V&x [4/^ - Az] } z sgn(0 



2 + 2z2sgn(e) 2 + 2z2sgn(e) 

ysgniO {3z {y^ - x"") + 2xV&^/W\ [l - ^'sgn(e)] } 



[1 + z2sgn(^)] 



V6v^ 



0. 



Despite the above complicated form, we can straight- 
forwardly see that using the explicit critical points pre- 
sented in Table HI the matrix Q acquires a simple form 
that allows for an easy calculation of its eigenvalues. The 
corresponding eigenvalues Vi (i = 1, 2, 3) for each critical 
point are presented in table IIIII 



Appendix B: Stability of the critical points at 
infinity 



Let us consider the Poincarc variables 



/9 cos sin?/;, Zr = psin^sini/', yr = /ocos'0, (Bl) 



Thus, by determining the sign of the real parts of these 
eigenvalues, we can classify the corresponding critical 
point. In particular, if all the eigenvalues of a critical 
point have negative real parts then the corresponding 
point is stable, if they all have positive real parts then 
it is unstable, and if they change sign then it is a saddle 
point. 



where p 



r = y/x'^ +y^ + z2, e [0,27r], and 



-f<V'<f|2il, [2^''. Thus, the points at "infinite" 
(r — > +00) are those having p ^ 1. Furthermore, the 
physical phase-space is given by the intersection of the 
regions 2{x^ + y^) < 1, + y1 — z^ > and the circle 



* We restrict the angle ip to this interval since the physical region 
is given by y > 0. 
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Cr. Point 


Exists for 


1^1 


1^2 


1^8 


A 


all 5,A 


3 
2 


4 •J V J yo£; ) 




B 


5 = 0, all A 


3 







C 


5 = 0, all A 


3 


3+y|A 





D 


5 = 0, A^ < 6 


A^ - 3 







E 


5 = 0, A^ > 3 









F 


A^ < 5 


-3 




-3-v'9-24V«"-A2C 


^ 2 


^ 2 


G 


A^ < 5 or 5 < 


-3 




-3-\/9+24V?^-A25 






J 


5 / and A = 


-3 


-3+V9-24e 
2 


-3-V9-24e 
2 



TABLE III: The eigenvalues of the matrix Q of the perturbation equations of the autonomous system (|24p . Points B-E exist 
only for 5 = 0, and for these points the variable z is zero and thus irrelevant. Although for these points the eigenvalue associated 
to the z-direction is zero, the stability conditions are obtained by analyzing the eigenvalues of the non-trivial 2x2 submatrix 
of Q. 



Cr. Point 


Exists for 




V2 


A± 


5 > 0, all A 


\ 


3 - \/965 - 24^/65 + 9 






3 + \/965 - 24^ + 9^ 




L± 


5 > 0, all A 


\ 


3 - \/965 + 24^65 + 9 






3 + \/965 + 24^/65 + 9^ 





TABLE IV: The eigenvalues of the matrix Q of the perturbation equations of the autonomous system (|B2|) . calculated for the 
four critical points at infinity. Since we are restricted to the invariant set j/r = 0, Q is a 2 x 2 matrix. 



- + 2/r + -^r — Ij that is it is the region limit p — > 1, the leading terms in (|B2[) are 

(a:„2/,,z,)e[-l,l]x[0,l]x[-l,l]|z2<2;2. + y2< 11 P 



a/I — p^d' ^ ^•^cos'0cot'0sin0, 
Performing the transformation (|B1[) . in terms of the 
Poincare variables the system ((M)) becomes: i/l — p^^' _^ J - Xcos 6 cos ip . (B3) 

x' = ^ ^ \ r I r r ^^j.'^ _ 2y'^ — l) X Since the equation for p decouples, it is adequate to inves- 

tigate the subsystem of the angular variables. Therefore, 
y/6£, [(2y^ - i) xf. - yl + l] Zr the singular points at infinity satisfy 

x'^ + yf. — 1 = ± cos 6, 

^J^Xyf- Zr = ±sin6', 

, 2y/6^yrZr {—x^ — — z'^ + 1) xl for some particular 9 e [0, 27r] that are determined as 

yr = y2 _ _|_ 2 (^2 -I- ?/2 _ j^^] foUows. Substituting ()B4p in ()B2|) and solving for 9' we 

obtain 



-yr [6x1 + 4:^/6£_ZrXr - 6yr + 3) 1 / ^ ^■ 

^ ^ ' 61' = --cos(26i) (^3cot6i + 2V6V4 



(B5) 



This equation admits the parametric solution 



- _ 2;2 + 1 _^ 

4^6? (x^+y,^- 1)^3 2V6?y,^z,^x. = ci + {4V6etanh-i(tan0) 
z2 + 2 (a;2 + y2 _ 2) a;2 _,_ y2 _ 2 

3 (2x4 - x2 - 2y4 + y2) +31n[cos(20)] - 6 In [3cos0 + 2y6esin( 



6^ {2xt + xr) + ' \, 7 — f—^^^.m 

Z (xf. + y^.-l) Thus, the singular solutions with 9 e [0, 27r] are those 

Transforming to spherical coordinates and taking the satisfying 9 = |,^, and cote* = -2w|$, which due 
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to (jB4[) lead to simple expressions for Xr,yr,Zr. These 
results are summarized in Table |lll 

In summary there are four physical critical points at 
infinity. The points K±, satisfying 9 = ^ and the points 
L± satisfying 9 = j. Thus, according to (|B1|) . these 

critical points correspond to the limit (j) ~^ ±cx), (j)/!! — > 
±oo satisfying the rate 

(ln0)' = VE- = V6cot9. 

z 

Concerning the stability analysis of the singular points 
K± and L±, we take advantage that they are located in 



the invariant set yr = and we examine their stability 
for the reduced 2D system Xr,Zr. The eigenvalues of the 
linearized 2D subsystem evaluated at K± and L± are 
displayed in Table IIVI while the results of the stability 
analysis are presented in Table [TTl 

Finally, we mention that mathematically there are two 
more critical points, namely the points P± satisfying 

cot 6' = — 2y^^. However, they have no physical mean- 
ing since they give a divergent VIde in ()26p . namely 
sign [8^ — 3] ■ oo, and thus we do not consider them in 
the cosmological discussion. 
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